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NOMENCLATDEE 


Cross-sectional area of stiffner 
Axial length of cylindrical shell 
Extensional rigidity 
Spacing of stringers 
Young’s modulus 
Resonant frequency (Hz) 

Objective function (wei^t of shell) 

Siear modulus 

Acceleration due to gravity 

Moment of inertia of stiffener about its centroid 
Moment of inertia of stiffener about middle 
surface of cylinder i^rall. 

Torsional constant 
Bending rigidity 
Spacing of rings 

Mass per unit area of cylindrical shell 
Number of axial/meridional half waves 
- Stress resultants in cylindrical shell 


^xy 

- Stress resultants in conical i^ell 
Qs 

- Number of circumferential full waves 

- Number of lings 

- Number of stringers 


Axial load 
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R 

R, 


R 


S 

S 

s 



u, V, w 
X, y, z 


Z 

z 


Mean radius of cylindrical shell 
Mean radius of cone at the .qma.~i i end 
Mean radius of cone at the large end 

S A 

Non-dimensional parameter, ' 

E t 1 
•w 

Rigidity moment 

E A 

Non dimensional parametei; g ' 

w 

Distance from cone ap^ along meridion to a 

point on shell middle surface 

Distance along meridion frcm cone apex to pm all 

base and large base of conical shell respectively 

Twisting rigidity 

Time (secs.) 

A 

Effective wall thickness of i^ell, — g 

Thickness of ring 
Ihickness of stringer 
Thickness of shell wall 

Displacements of a point on shell middle surface 
in X, y, z or ©, s, z directions, respectively 
Displacement amplitudes 

Orthogonal coordinates as defined in Eig* 2.1 
(x and y lie in the middle surface of cylinder) 
Curvature parameter, — (1 - 

Distance of centroid of stiffener to middle 
surface of shell, (positive if stiffener 
lies on external surface of i^iell) 
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©, s, z - Orthogonal coordinates as defined in Fig. 3-1 
(Q and s lie in the middle surface of cone) 


Greek Symbols 


r 

& €: €. i 

^9, s, 9s 


Wave length pai^eters or angles (Fig. 3*1) 
%)ecific weight 
Middle surface strains 




rr 

X 

/ 


00 


- Poisson’s ratio 

- Potential energy or change, in 

potential energy 

- Curvature 

- Mass density 

- Circular frequency (rad/ sec) 


Subscripts 

c - Cylinder or cone 

r - Bing 

s - Stringer 

CO - Inertia load 

w - Shell wall 

A subscript preceded by a ccanma indicates 
partial differentiation with respect to the 
subscript • 



SINOPSIS 


In the present woric, the design optimization of 
axially loaded, sii!5)ly supported stiffened cylindrical and 
conical shells for minimum weight is carried. The design 
variables are thickness of shell wall, thicknesses and 
depths of rings and stringers, number/ spacing of rings and 
number/ spacing of stringers. Natural frequency, overall 
buckling strength and direct stress constraints are con- 
sidered in both cylindrical and conical shell design 
problems. In the case of cylindrical ^ell the local 
buckling constraints are also included. 

The design of cylindrical shell is carried with 
three different ccmbinations of stiffeners. In the case 
of conical shell, an expression for the buckling load is 
derived and optimum designs are obtained with rectangiiLar 
stiffeners. For all the examples, the independent effects 
of behaviour constraints are also studied separately. 

The optimum designs are achieved with one of the standard 
nonlinear constrained optimization techniques (Davidon - 
Fletcher - Powell method with interior penalty function 
foimiiLation) and few optimal solutions are checked for 
the satisfaction of Kuhn-Tucker conditions. 



CHAPTER I 


INTEODUCTION 

1.1 IMPORTANCE OF THE PROBLEM: 

One of the most cominon structural eloments used in 
modem airplane, missile, booster, and other aerospace 
vehicles is the thin conical or cylindidcal shell. In 
aerospace vehicles these shells are used as Interstage 
adapters, which mate ttie space-craft to the rocket vehicle. 
Biese shells are weak in resisting the frequently encoun- 
tered high compressive stresses. Hence to achieve higher 
efficiencies, it is necessary to design the shells with 
high buckling strength-to-weight latio while satisfying the 
minimum natural frequency and yield strength requirements. 
One of the common means of achieving the above requirements 
is to stiffen the shell by longitudinal (stringer) and 
circumferential (ring) stiffeners. Thus the design probloE 
becanes one of finding the skin thickness of the ^ell, 
sizes of the stiffeners, and the number of rings and strin- 
gers. 


1 .2 STATE OF ART 

Overall buckling and free vibration analyses of 
stiffened cylindrical and conical shells have been of 
interest to designers for a long time. Because of the 
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canplexity in the analyses, earlier designs are "based on experi- 
mental data and local huclding equations* But the data has 
vride scatter and hence one has to make a conservative esti- 
mate of the critical values of buckling load and natural 
frequencies. The designs obtained from these critical 
values may not be optimum. However, some woik has been done 
by using approximate analyses. 

The buckling analysis of cylindrical shell has been 
made by many investigators [1 - 12 ]. Serpicof 1] and 
V/eingarten et al. [ 2, 3) ^ ] have studied cylindrical and 
conical shells without stiffeners. Later tiie effect of 
stiffeners has been considered in the analysis by Block 
et al. t5l and Singer et al. [6]. In Ref. [51 and [6] , 
the effect of stiffeners has been averaged over the shell 
surface. A few research workers [8, 9? 12] have accounted 
for the discreteness of stiffeners. 

Free vibration analysis of stiffened cylindrical 
[ 13 - l5] and conical [16, 17] shells has been carried by 
few research workers. In Ref. [13] an expression for the 
natural frequencies has been derived, considering the 
effect of stiffeners by smearing technique* ■-Late2‘'Efle 
and Sewall [ l4j have analysed cylindrical ^ells account- 
ing for the discreteness of the stiffeners. Crenwelge 
and Muster [ 17 ] , using the energy approach, indicated a 
method to determine the natural frequencies of simply 
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supported conical she3.1s in unstiffened, ring- stiffened 
or ring-and-stringer stiffened configurations. 

A wealth of literature [ 18 - 38] has been devoted 
for the optimum design of cylindrical shells subjected to 
general and local buclLling constraints. Before 1967} the 
studies 1 21 - were based on the assumption of simulta- 
neous failure modes and the optimization was achieved by 
parametric studies, ilfterwards optimum design has been 
achieved by the use of constrained optimization techniques. 
Pappas and Amba Eao[30 lhave considered the designs with 
conventional as well as spiral stiffeners. In 1975} 

Of 

Fronowicki et al. [38 ] presented the design shell with T - 
ring stiffeners taking three different objectives subjected 
to vibration constraint. The three objectives are (i) mini- 
mum weight of the stiffened shell, (ii) maximum separation 
of the lowest two natural frequencies, and (iii) maximum 
separation of 'the lowest two natural frequencies which have 
primarily axial contents. 

The design of stiffened conical shells has been 
attempted by only few research woiiiers. In 1966, Bums 
[39 }has achieved optimum design of ring- stiff ened shells 
by parametric study. The loading considered was hydro- 
static compression. Later Ihomton[401 and Herald et al. 
[41 ] have used nonlinear piogramming techniques for the 
design of stiffened conical shells v/ith buckling and 
yield stress constraints. 





1.3 OBJECT MD SCOPE OF PEESENT IJOEK 

The object of this woik is to present a iinified 
approach for the design optimization of simply supported 
stiffened cylindrical and conical shells for minimum weight. 
*hxial loading is asstimed for both the sliells. For conical 
shells, an expression for ciitical buckling load is deriired. 
In the case of cylindrical shell, examples with tiiree 
different combinations of stiffeners are considered. An 
example -idiich was solved in Ref. f 26 j is also solved with 
the constraint set derived in this work and the results are 
compared. Only one t3^e of stiffeners (rectangular) are 
considered in the case of conical shell.. In all the exam- 
ples, the effects of frequency and buckling constraints are 
studied separately and the results are compared, Eie 
examples are formulated as standard constrained optimiza- 
tion problems and are solved by gradient method (Davidon - 
Fletcher - Powell method) with interior penalty function 
formulation. Finite difference (backward differences) 
technique is used for the gradient calculations. Tire opti- 
mum solutions have been verified for the satisfaction of 
Kuhn - Tucker conditions in few cases. 
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1.4 ORGANIZilTIOH OF THE THESIS 

The thesis is divided into six chapters. In 
second and third chapters, the buckling and vibration 
analyses of stiffened cylindrical and conical shells are 
presented respectively. The formulation of the optimiza- 
tion problem and solution procedure are given in chapter 
four. Results are presented and discussed in chapter five. 
Chapter six gives the conclusions of the present study and 
recommendations for future work. 



CHAPTER 2 


• ANALYSIS OF CILINDEICAL S-IELL 

Bie buckling and vibration analysis of stiffened 
cylindrical ^ell, presented in tiiis chapter, lias been 
originally given by Block et al.[ ^1 and Mikulas et al.t 13 I 
based on energy principles. Their a.nalysis is based on 
small deflection theory and includes eccentricity Cone - 
sided) effects of the stiffeners. The results of this 
analysis have been shown to be in good agreement \.ith other 
theoiiesand experimental results. The stiffened cylinder 
is considered to be composed of an isot]?opic shell stiffened 
by uniform, equally spaced rings and stringers, all having 
elastic material properties. 

The following assumptions have been made in this 

analysis 5 

(i) the rings and the stringers are spaced closely 
so that their elastic and inertial properties 
can be averaged over the stiffener spacing, 

(ii) the usual Donnell type assumptions are used to 
specify the displacements in the shell, idiere 
as the stiffeners are treated as beam elements 
with stiffener twisting accounted for in an 
approximate manner. 


\ 
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(iii) in cases where rings and stringers lie on the 
same surface of the shell, the effect of the 
joints of the stiffener frame is ignored. 


2.1 BUCICLING ANALYSIS 

For the coordinate system diown in Fig, 2.1, the 
strain- displacement relations in the shell due to hucltling 
displacements can he written as 


with 


e 

X 


e 


xy 



- 2 W 


»XX 


” ^ ^’yy 

- 2z w 


xy 


( 2 . 1 ) 


ey - "^Jy R 
®xy = ^»y + ^»x 

where u, v, and w are the displacements induced hy 
hucld-ing . 


If the stiffeners are assumed to behave as bean 
elements, the stiffener strain^- displacement relations 
can be written as 


e = € — 2 w. 

X X »xx 

s 

7j. = - Z "'^Jyy 


( 2 . 2 ) 
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vjhere the subscripts s and r are used to denote strin- 
gers and rings, respectively, iScuation (2,2) specifies 
that buckling displacements are such tlmt the strain varies 
linearly across the depth of the stiffner, yet satisfies 
the compatibility of displacements between the stiffener 
and the surface of the shell to whidi it is fastened. This 
compatibility requirement is the source of eccentricity or 
one-sided effect in the buckling of stiffened shells. 


The strain energy of the isotropic shell ^ can 
be expressed in terms of middle surface displacements as 


2irR 


2 (1 


I tuL + 



The strain energies of stringers u g and of 
sings corresponding to the buckling displacements 
can be expressed as 



a 


9 


TT 

S 


1 

2 d 


2 ttR 



/ 





2 z„ E A u, w. 
s s s »x >xx 


+ lo 
s s 


'^’xx ^^^xy J ^ 


(2.4) 


and 


2 X R a 


1 

2 1 

/ 

o 

c 


+ -) 

R'' 

''>yy ''?yy ^ 


J-i- 

r r r 


(2.5) 


vAiere d is the stringer spacing and 1 is the ring 
spacing. The subscripts s and r are used to denote 
stringer and ring properties, respectirely. z denotes the 
distance from the centroid of a stiffener to the middle 
surface of the shell as shown in Fig, 2.1, and is positive 
if the stiffeners are on the external surface of the cylinder, 
lo is the moment of inertia of a stiffener about the middle 
surface of the shell. 


In Eqs. (2.4) and (2.5), the energy contributions 
of the stringers and the rings are averaged over the cir- 
cumference and the length of the cylinder respectively''^ 
and hence, discreteness of the stiffeners is ignored. 

A rigorous derivation of the buckling equations 
from the potential energy of -the loaded shell reauires the 
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use of nonlinear theory and techniques. However, when the 
effects of prehuckling deformations are not considered 
explicitly, the potential energy of the external forces 
causing huddling can be derived as follows. 

As the middle surface of the shell bends due to 
buckling, no membrane stresses are produced but the edges 
get closer to each other ^iiich give rise to inplane dis- 
placements. The inplane displacement of an element of 
length dx in x - direction can be obtained from Fig. 2.2, 
as 

Au = ds-dx = '5 ^ (2. 6a) 

Consequently, the work done due to the force Iw 
acting on an element of width dy of shell element can be 
expressed as 

A u = I ( If ) dx (2.6b) 

Similarly, considering the work done by the forces 

and N , the potential energy of the external forces is 
y 

2 IT R a 

’'l = - i + 2 

o o 

+ n w? ] dx dy (2.6) 

where N is a stress resultant (positive in compression) 
obtained by considering the shell and stringers to be 
loaded with a uniform normal stress in tire x - direction. 




FIG.2-2 STRAIN, £x^DUE TO LARGE 
DEFLECTION (w; 
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Similarly is a stress resultant (positive in compression) 
obtained by considering the shell and rings to be loaded 
xith a normal stress in the y - direction, 

Ihe total potential energy of the loaded stiffened 
cylinder, ^ can be written as 

”1 = ’c * ''s + ’’r ”1 <2.7) 

The application of the principle of miniinum poten- 
tial energy ( 6 ir ^ = 0) yields the equilibrium equations 
and the boundary conditions corresponding to arbitrary 
variations in u, v and w respectively. The displacement 
functions which satisfy the simply supported boundary con- 
ditions [v (x=0,y) = V (x=a.,y) = 0 and w ,Cx=G,y) = 


V (a:z£L,y) 

u = 

= q 

u 

] are 
cos 

mTT X 

a 

cos 

m 

R 


Y = 

Y 

sin 

z 

a 

sin 

ny 

H 

(2.8) 

w = 

W 

sin 

m TTX 

a 

cos 

SZ 

R 



where m denotes the number of longitudinal half v/aves 
and n indicates the number of circumferential full xravss 
in the buchled pattern. 

If Eq. (2.8) is substituted into the equilibrium 
equations obtained from = 0 wiUi = 0, the 

existence of nontrivial buchling displacements require 
that the determinant of coefficient matrix of u, v and w 
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vanish. Biis condition reduces to the stability 


eoiiation : 




/-I • a , 2 "^s 2 ^ ^ 

-tB } + in -- 3 - rr- + m B 


=: m 


I. 


^ d D '*' 1 D ■ ) 21 B 


-1 2 z" 

^ ^ nr I 


1 + S A + R A + s R A 

— £__ rs 


(2.9) 


■where 


= 1 + 2 (1 - S^M ) ( -^ ) (1 + 6^)2 (h) 


Z_ 2 


1 + 2 a 2 ( 6^ - y) ( -| ) + (1 +6 2)2 ^ ^ 


1-v2,2„2b2(1-„2) 


^ r -I It 2 . _ -2 / . s , ^"Pv^ 


+a 6 [1 -y + 2 (1 +y) ] ( ^) 


20.*^ 6*^ (1 ^y)2 ( h , ( |s , j2 


[ 2 (1 + n) + (-I . p2) j 


(1 + B ) + 2 B^ (1 +y ) ( E + S ) 


+ (l-y^)[S + 2B^ES(1+y)+B^E ] 


with 


R^ t^ 
w 


D = 


“12 (1 - y2) 



1 ^ 



-jj _ R _ n a 

^ ~ m irH 

In Eq, (2.9), the specified comhination of axial 
and circianf erential loading must be minimized numerically 
for integral value of m and n in order to compute the 
critical buclfLing load. 


2.2 FREE VIBEA.TION ANALYSIS 

In the case of free vibration analysis, the total 
potential energy of the stiffened cylindrical shell 
is composed of strain energy of the shell, -jr^J strain 
energy of the stringers, strain energy of the rings, 

T^, and potential energy of the inertia loading, ir^ , 

The expressions for ir g and it p are the same as 

defined in Section 2.1 ( Eqs. (2.3), (2.4) and (2.5) ]. An 
expression for can be ob'bained as follows. 

If the shell is undergoing simple harmonic motion 
with circular frequency m (Inplane inertia neglected), and 
w (x, y) is the deflection shape at the time of maximmi 
deflection, the potential energy due to Inertia load can 
be written as 

2ir R a 

IT = - ^ f / M 0)^ w^ dx dy (2.10) 
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where M = , is the 

average smeared-out mass per unit area of the shell. The 
quantities P^, and are the mass densities of 
cylinder wall, stringers and rings, respectively. 

The application of Haiailton* s piinciple or 
principle of minimum potential energy in equivalent sfc’tic 
case ( ^ = 0) yields three equilibri'um equations. On 

substituting the displacement solution, given by So. (2.8), 
in the equilibrium equations, a set of homogeneous linee^r 
equations is obtained in u, v and w. Applying the con- 
dition for nontrivial solution, an expression for frequency'' 
is obtained as 


A F k p 2 V E I 




d D 


+ B ^ ( 


G J G 1 K 

s s . r r \ A 


d D 


1 D 


) +6 


E_ I_ 
.JL^ 1 
1 D 


^ 12 Z‘ 

+ -IT- 
ir 


1 + SA +RA +SE A 
[ ® 


( 2 . 11 ) 


where 


F is the frequency in cycles/ sec, 

o ^z 2 o2 

Ag = 1 + 2a 2 ( -^ ) ( 0 2 A ^ ) d+B"^) 

Ay = 1 + 2 n^ ( ^ ) (1 - ^) + n^ ( -^ ) (1 + 6^) 
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^rs - 1 - + r? ci^ (1 - P^) + 2 (1 + |i)i ( -^ )^ 

h- p 5 z 2 

+ n [ 1 - !.■" + 2 6“^ (1 + p) , { ^ ) 

+ 2 (1 - p2) [ ( !| ) + ( fj ) , 

+ 2 n“^ (1 + p)2 ( fs ) ( !| ) 

and 

2 — 

•A, Z , S, R, a J and B are same as defined in Sq, (2.9). 

In order to get the minimnn frequency, the right 
hand side of Eq, (2.11) must be minimized numerically for 
integral values of m and n. 


2.3 STRESS MALYSIS 

(i) Stress in Shell Wall: 

For an axially loaded stiffened cylindrical sliell, 
the stress induced in the shell v/all can be tahen apqroxi— 
mately as 


a 

w 


P 

2 TT R t 


( 2 , 12 ) 


'Where P is the axial load, 

R is the mean radius of the; shell, 

and t is the effective thickness of the stiffened shell 

_ A 

in axial direction f43l, (t = — f+t), 

^ d W 
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(ii) Local Buckling Stresses s 

If the stringer spacing is anall relative to the 
radius of the cylinder, classical flat plate theory can be 
used to predict local instability in the cylinder wall 
(panel buckling )[ 22]. 

The critical buckling stress in the cylinder i/all 
is given by 


"fc 2 

0 cr^ = t 4 Eyi2 (1 - y^)] ( ) (2.13) 

Similarly, local instability in rectangular 
stringers can be predicted using the following equation 

[22] J 

(Tcr^ = [ 0.5 ir^ E^12 (1 - y^) ] ( ) (2.l4) 

where 

t 2 is thickness of stringer and 
1 is spacing of rings. 


If stringers are considered to fail in the column 
buckling mode, the critical load is given by 


P 


cr 


0.5 




( 2 . 15 ) 


where 1 is length of the stringer between two rings 
(equal to ring spacing), 

and I is moment of inertia of stringer about its 
s 


cent roi dal axis 



CEIAPTER 3 


ANALYSIS OF CONICAL fflELL 

In this chapter, the buckling and free vibration 
analysis of stiffened conical shell is presented. The shell 
is assumed to be loaded with an intensity of P (load per 
■unit length of circumference) at the snail end in the axial 
direction. For this case, an expression for the buckling 
load intensity is derived in Section 3»1» In Section 3»2, 
the free vibration analysis is presented as per the method 
indicated by Crenwelge and Muster 07 1 • The results of 
free vibration analysis have been shown to be in good agree- 
ment with other theoiles and experimental results [l7j» 

The following assumptions are made in the buck- 
ling and vibration analysis of conical shells 

(i) shell displacements and displacement gradients 
are small so that linear shell theory is 
applicable, 

(ii) rings and stringers are closely spaced so that 
their elastic and inertial properties can be 
averaged over their respective spacings, 

(iii) the stiffeners are thin, the interactions at 
the ring and stringer joints may be ignored, 

(iv) stiffeners are twisted about their centroids. 



3 . 1 BUCICLING ANALYSIS 
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Consider a tliln truncated conical shell, whose 
geometry and coordinates are shown in Hg, 3,1. Bie 
strains at a distance z from the middle surface of the 
shell are given by 


€ =: £ — Z Y 

S S ^ S 


£0 = Eg - z X0 

0 ~ ^ S0" ^ s 0 

For the rings 

*' z Y„ 

^ c 0 ^^0 


S0 


and for the stringers 


(3.1) 


(3.2) 


e - z X 

s s 


w 


»S0 


2 w, 0 


u, tan a 


s0. 


s cos„ 2 

“ s cos a 


(3.3) 

2 u tana 


■idiere 


e = V- 
^s ^ s 


e« = 


'0 
^se 
5(s 


(u, ^ cos a + V + w tan a) / s 
u,g - u/s + V, 0 / (s cos a) 
w, 


ss 


w 




^ ^ ^^*06 " % ) / (s^ cos^ a) 
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30 S COS 


a „2 


s cos a 


u, tan a 

—IS u tan a 

2 s 2 J 

s 


and u, T and w indicate buckling displacements. 

Tbe strain energies of cone, stringers and rings 
in terns of strains are given by 


o 


S }a.V fv2,„2 


^ ^ ^0 X s + 2 li ^ X s X 0 + 2 C1 - y } 1 


s cos a d 6 ds 


(3.^) 


Sp 2 ir 

2 ^ ^ f ^s^s^ ” 2 „ + K^X 


o 


S^sXs 


^ s e J ^ cos a de ds 
s 

Sp 2Tr 

? J /tD^^6^-2S. ■ 


( 3 . 5 ) 


6 - 2 Sj, S3X0 + KJ.X, + Tj.x^g 1 


S.^ o 


where 


s cos a de ds 


E t / (1-y 
w w w 

S A / b 
r r r 

A / b 
s s s 


( 3 . 6 ) 


■^r ^r ^ 
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E A z / b 
s s s s 

E„t3/ 02 (1 -u/), 

\ «r ^ / t-r 


= E^ (I^ + z/) / b„ 

s s s ss ' s 


T = G J / b 
r r r ^ r 


T = G J / b 
s s s s 


and b =3$ cos a 
s 


Tlie potential energy of the external forces 
causing buckling is given by 

T 

\ = - 2 t®9''?e^2Hg3W,gW,3 +K^<3l 

O 


s cosct d6 ds 


(3.7) 


For an axial load (compression) of intensity P 
(load per unit length around the circimference), the stress 
resultants and N_ are given byt ^2 ] 

0 ^ So 3 


(3.8) 


s sin c( 


The total potential energy of the shell can be 


expressed as 

V _ iTj-TT + ir + Tr 

T - c s r 1 


(3.9) 
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‘Ihe simply supported end conditions are 
u (s^, e) := u (s^, © ) = 0 

and w (s^, e) = v (s^, P) = 0 

The displacements of the middle surface of the 
shell, which satisfy the boundary conditions [ Sq. (3.10) ], 
are given by 


11 

— 

2 

sin 

=: U 

s 

Y 

— ^ 

2 


= Y 

s 

cos 



2 

sin 


= W 

s 


where 
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(3.11) 


m is the number of meridional half-waves, 
and n is the number of circumferential full waves. 

2 

The term 3 is proportional to the cross sectional 
area enclosed by the shell at a distance s fran the cone 
apex. It is used to shift the point of maximim displace- 
ment towards the large end of the cone, a phenomenon that 
occurs in reality due to decrease in stiffness of the 
shell with increasing radius [ 17J . 


Substituting the solution [ Iq. (3*11)1 iri 
Sq, ( 3 . 9 ) and applying the condition for minimum potential 
energy, a matrix equation is obtained as 
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where are defined, in Appendix A, 

a,nd is the critical load intensity (load per unit 

length along the circumference) given by 
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3.2 FREE VIBRATION ANALYSIS 

Ihe total energy ttj of the stiffened conical 
shell, ■^^h.ich is in a djmainic state, is the sun of stirain 
energies of the shell stringers rings ir^, and 

potential energy due to inertial_ load (kinetic enei^y) 

^v-"' ^c’ ^s "’’^r same as defined by Eqs, 

(3.^), (3.5) and (3.6) respectively, n is given by 
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The displacement solution for this analysis is 

assumed as 

u = u S^ sin (s - s.^) cos n0 

V = V e^“^ S^ cos \ (s - s^) sin nS (3.16) 

w = w e*' S sin (s - s.j) sin ne 


where = niir / (s^ - s^) , 

m is the number of meridional half-\-javes, 
and n is the number of circumferential full waves. 


Substituting the solution (Eq. 3.16) in the total 
energy expression (Eq, 3*15) > and applying the condition 
for minimum potential enei^, the equations of motion for 
the vibrating shell are obtained as 
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where are given in Appendix A, 

'Ihe characteristic equation corresponding to 
Eq. (3.17) is 
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(3.18) 
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where 
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Of the three natural frequencies obtained fr<xi 
Eq. (3.18) for each combination of m and n, only the 
lowest (v/hich is associated vjith the ladial motion of the 
shell) is of interest. 


3.3 STRESS MALYSIS 

For a given axial load (compressive) of intensity 

P, the resultant stress [Er, (3.8) ] N„ is inversely 

s 

proportional to s. Hence ‘bhe maximum occurs at the 
small end (s = s^) of the shell and the magnitude is 

given by 
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CHAPTER h 


FOMDLATION AND SOLUTION OP OliS PBDBLM4 

This chapter deals with the formiiLation of optinrum 
design probleiiis as standard nonlinear prograimning problems 
and also the techniques used for solving these problems# 

The design problem of stiffened cylindrical shell is for- 
mulated in terms of objective and constraint functions in 
Section ^*1, Section 4.2 deals with the formulation of the 
problem of design optimization of stiffened conical shell 
under buckling, natural frequency and direct stress cons- 
traints. The procedure forjsblving the design problems 
is described in Section 4,3# 

4.1 FOR^IULATION OP CYLINDRICAL SHELL PRDBLai 
4.1.1 Ob j ective Function 

In the design optimization of ejCLindrical shell, 
the objective function normally considered is weight. For 
the present problem, the sum of the individual weight of 
shell wall rings, Np and stringers V/^ foims the objec- 
tive function. Thus the total weight, f, of the siiell is 
expressed as 

f = afTE a + a O.la) 

where R = mean radius of cylinder, 
t = thickness of shell wall, 
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a = length of shell, 

/ , / _ = specific weights of shell -trail, rings and 

w X s 

stringer materials, respectively, 

A^, Ag = cross sectional areas of ring and stringer, 
respectively, 

n , n = n-umber of rings and stringers, respectively. 
X s 


^.1.2 Design Vaidables 

In majority of the applications, ^ells are used as 
components of a system. For a shell to be compatible ^jitii 
the mating parts, overall dimensions such as mean radius and 
length of the shell are fixed beforehand. Consequently, 
the parameters that are taken up as variables for the optf-Kum 
fih el 1 design are shell -vra.!! thickness, cross sections of 
stiffeners, number/ spacing of rings and number/ spacing of 
stringers. In the present problem, spacings of si-iffeners 
are treated as continuous variables. The objective function 
Eg. C^.1a) can be rewritten in terms of spacings of 
stiffeners as 

f = R \ a + 2nR + Ag a ^ i 

( 4 -. 1 ) 

where the design vsiriables are given by 


t = "tiiickness of ^ell -wall 
•w 

= thickness of ring 
b^ = depth of ring 
t^ = thickness of stringer 



= depth of stringer 29 

1 = spacing of rings 

d = spacing of stringers 

>+. 1.3 Behaviour C onst raint s 

Shells, as they are used in space-craft and aero- 
space industry, have to satisfy certain behaviour constraints 
like natural frequency, overall buckling strength, local 
buckling strength and direct stress. These constraints are 
formulated as follows. 

(i) Constraint on Natural Frequency 

This is one of the important constraints which must 
be satisfied to avoid failure by resonance. liesonancc 
occurs when disturbing frequency is equal to one of the 
natural frequencies of the shell. Hence minimum natural 
frequency must be constrained to be not less than tiie ciis- 

turbing ones. 

Theoretically, the shell has infinite natural 
frequencies, which are given by Eq, (2.11) for different 
integer values of m greater than or equal to one and n 
greater than or equal to zero. In practice the minimum of 
these, which occurs for finite values of m and n, are 
constrained, is the combination of m and n which 
gives minimum frequency is not known befor^and, freauen- 
cies have to be computed for all combinations of m and n 
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until minimtim is obtained. For verification, finite upper 
bounds for m and n are assumed and natural frequencies 
are computed for a number of ^ells with different stiffe- 
ners. From these results, it is observed that minimum 
frequency in each case occurs al\i?ays at m equal to one 
only. Hence m can be taken as one and n has to be 
varied to find the minimum frequency. In an optimization 
routine, this minimum has to be found a number of times, 
which requires a lot of computational time. Hence the 
following procedure is adopted, 

A design vector of independent variables is assumed 
and the frequency spectrum is obtained at the initial o.esign 
configuration of the shell. In this spectrum let the value 
of n corresponding to minimum frequency be n’ . -*^s the 
design vector changes in optimization routine, the value n 
which gives minimum frequency will also change but it is 
expected to be in the neighbouihood of liie starting value 
n* . Hence frequencies at the values of n equal to n’ - 1, 
n' and n’ + 1 are constrained as 


F. ^ F . 

1 ' min 


i — 1 , 2, 3 


(^-. 2 ) 


where F^^^ are taken as F^ = F (m - 1, n = n‘ - 1), 

F = F (m = 1 , n = n') and F. = F (m = 1, n = n» + 1), 
2 

and ^min lower bound on frequency. 
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(ii) Constraint on BucltLing Strength 

In practice, ^ells are subjected to different 
combinations of axial, external and internal loads. These 
loads sometimes lead to the failure of the shell due to 
bucltiing before the induced stress reaches the yield 
stress. Hence, a constraint on the buckling strength is 
necessary to overcome this failure. 

A procedure similar to the case of natural fre- 
quencies is adopted in forming iiie constraints on the 
budding strength, except that m is not kept constant in 
this case. If the values of m and n corresponding to 
the minimum buckling load are m” and n” at the starting 
design vector, buckling loads at different combinations of 
m' (viz., m” - 1, m” and m" +1) and n (viz., n” - 1, n 
and n” + 1 ) are constrained as 

Id 2 /TR t 2:P applied X Factor of safety, t 

i = 1, 2, ... 9 C^*3) 

>±Lere N are given by Sq. (2.9) with Id = 0 

X. J 

_ ^ A 

and t =: t + . 

W U. 


(iii) Constraint on Local Buckling Strengths 


Besides the constraints on frequency and ovemll 
buckling strength, it is safer to have constraints on 
local buckling strength also. This need arises due^ to , 


Oc.M I i % ’ V; , ‘ 

has* ^ 


i5284- 
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local imperfections and inaccuracies in fabrication. 
Jailure may occur because of local buckling of the shell 
wall or buckling of stiffeners. 

considering the panel between any two stringers, 
the expression for critical buclriLing stress is given by 
Eq. ( 2 . 13 ). Ihe constraint on this critical stress is 
expressed as 


^ ^^applied ^ ^a<^tor of safety/2/TB t) (4.4a) 

Similarly the constraint on critical bucltling 
stress in rectangular stringers given by Eq, (2.14), can 
be expressed as 


CTcr^ ^\pplied ^ factor of safety/2 /IE t ) (4.4b) 

In addition to the above, the stringer also may 
fail due to column buckling. The constraint on the buck- 
ling load Eq, (2.1^) can be expressed as 


cr 


^ (P 


applied • 


Factor of safety) 

(4.4c) 


(iv) . Const laint. bn Yield Stress 

The induced stress in the shell wall is given 
by Eq. (2.12) sliould be less than the allo^jable stress 
of the material, i,e., 

^ ^applied/ 2 ^ ^ ^ ^ (Yield stress/Factor of safety) 

(4.5) 
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Besides the above constraints, each design 
Variable is constrained by a lower bound. 

4. 2 FORMULATION OF CONICAL SHELL PROBLEM 

4.2.1 Objective Function 

In the design of conical shell also, the total 
weight of the shell is talcen as the objective function for 
minimization. The total weight, f, wliich is the stun of 
weights of skin (¥^), rings (W^) and stringers (¥^) is 
given by 

f = + ¥2 + ¥^ (4.6) 

where 

Wj = at [ (aj_ + (1 - 0.5)bj. sln(90-(i)i]A^Yr 

and Wj = (Sj - S3_) 

4. 2. 2 Design Variables 

Shell thickness, cross sections of stiffeners, 
number of rings and number of stringers are considered as 
variables for the optimum design of conical shell. Similar 
to the case of cylindrical shell, the radii of the conical 
shell and the angle a (Fig. 3.1) are assumed to be fixed 



beforehand. 

The design variables are 

^w = 

thickness of sliell wall 

ti 

thickness of ring 

b^ 

depth of ring 

t^ 

thickness of stringer 

b2 = 

depth of stringer 

^r = 

number of lings 

^s = 

1 ^ ^ 

number of stringers. 


^.2.3 Behaviour Constraints 

In the case of conical shell constraints on freauenc 
Buckling and yield stress are imposed as in the case of cylin- 
drical shell. These constraints can be stated as follows. 

(i) Constraint on Natural Frequency 

Unlike in the case of cylindrical shell, inplane 
inertias are considered in the vibration analysis of conical 
shell . Hence for a given set of values m and n, there 
exists three natural frequencies. Of these three, the lowest 
one wliich corresponds to the radial motion of tiie shell is 
of interest for any set of values m and n . The minimun 
of such lowest ones sho'old not be less than a p re dete mined 
constant. To reduce the computational time required in 
finding -fche minimum freq,uency, the following procedure is 
adopted. 
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It is observed that minimim frequency occurs for 
m equsl to one, and hence m is set equo.1 to one. The 
value of n which traps the miniiaum is sunposed to be 
finite. This finite value is taken as six wliicli is the 
case with the starting design vector for the example prob- 
lem considered. The resiILting constraints are expressed as 



^ F . 

• min 


where F^^ = F (m = 1 , 

defined by Eq. (3»18), 


i = 1, 2, ... 7 (^.7) 

n=i-1), i = lto7 F is 


(ii) Constraint on Buckling Strength 

Conical shells, subjected to axial loads and or / 
ciremf erential loads or both, may fail due to lack of 
buckling strength. Hence, buckling strength st.ould be more 
than a predeteimined value. This constraint is foimulated 
as described below. 

Buckling load, as expressed by Sq. (3.13), is a 
function of the axial and circumferential mode numbers m 
and n. Hence, a procedure similar to ohe one described for 
Eo.. (^.3) is adopted, to express this constraint as: 


P ■Sfc P 
cr_. ^ applied 


where 


CPi 


i = 1, 2, ... 9 

are defined by Eo. (3*13) 


(V.8) 


applied- 


is the applied load intensity, 


and 
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(iii) Constraint on Yield Stress 

To avoid failure of the material of shell vrall hy 
yielding j the induced stress in the shell given hy Bq. (3*19) 
should be less than the allo\7able stress, thus 

^apnlied ^ Tiela stress 

-—Hi ^ ^ 

t sino? Factor of safety 

In addition .to the above behaviour constraints, 
each design variable is constrained by a lower bound value. 

4-. 3 SOLUTION PEDCEDUEE 

The design optimization problems formulated in 
Sections 4,1 and 4,2 can be stated as standard nonlinea,r 
programming problems as follows s 

Find X which minimizes the objective function 
f ( X ), subject to the constiaints 

g. ( X ) S' 0 j = 1, 2, p (4.10) 

J 

where X is the vector of design variables. 

These problems are solved by gradient method 
(Davidon - Fltecher - Powell method) with interior penaJLty 
function formulation. For one dimensional minimization, 
cubic interpolation method is used [44] . As the function 
to be minimized is not in a simple and explicit form, 
finite differences (backward differences) technique is 
used for calculation of the gradient. Since the 
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Kuhn - Tucker conditions are necessarily to be satisfied at 
the optimum point, the Lagrange multipliers ( ^ ) are 
obtained from the gradients of active constraints (g., 
3=1, ... p') and the objective function (f ) [ 4 if]aB 


A = 

{[G] 

> 

J 

where 

=1 

Sx 

^S2 

Sg 

[g] = 

ag2 
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• 
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asi 
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d Sp 
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A = ^2 > 

i^P'. 
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7 - i Ai. 

-C = < “ 
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I df 

■ ■ ■^9 . 


and q is number of variables. 


The condition for any design vector to be a 
local optimum, the Lagrange multipliers computed at iiiat 
vector Eq. (i+.ll) should be positive. 



CHAPTER 5 


KESTJLTS AND DISCUSSION 

The optimization problems forrauLated in chapter 
four are solved and the nimierical results are presented in 
tills chapter. 


5.1 CILINDEICAL SHELL 

With a view of studying the effects of different 
sets of constraints and also the influence of type of 
stiffeners on the optimum solution, several cylindrical 
shell optimization problems are solved. The different 
problems considered are summarized in Table 5»1* 

TABLE 5.1 


TvnpR of stiffeners used fi 

ProIxLems solved 

for 

stringers 

0 for 

5 rings 

0 

0 

0 

0 with fre- 5 with oveiMJ winn all 

0 quency and 0 all buck- 5 cons— 

5 side cons 0 ling and 5 traints 

1 traints 0 side cons4i 

6 2 traints 5 

Rectangular 

Rectangular 

y 

y 

y 

Z- section 

Tubular 

y 

y 

y 

Z- section 

I- section 

y 

y 

y 


The data used in all these problems is given 
in Table 5.2. The numerical results are shown in Table 
5.3 to 5.5. The computer time taken for the solution of 
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an optimization problem ranged from ^..ll minutes on 
IBM 704I|- computer, 

5,1,1 Effect of Constraints 

For studjring the influence of different sets of 
constraints on the problem, the same starting point cuid 
side constraints are used for any particular combination 
of stiffeners. 

I'he following observations can be made from the 
numerical results obtained s 

(a) VJhen the behavioural constraint set consists of only 
natural frequency restriction, the frequency cons- 
traint has become active for all combinations of 
stiffeners. Further in all these cases, the spacing s 
of rings and stringers approached their respective 
upper bounds, while the thicknesses of skin, ring and 
stringer and depth of stringer have reached their 
respective lower bounds at the final solution. 

(b) When overall buckling is considered as the only beha- 
viour constraint, it has been found to be critical 

at the optimum point for any set of stiffeners. 

Here also the spacings of rings and stringers approa- 
ched their upper bound values and the shell thickness 
reached its lower limit. For all types of stiffeners 
all the side constraints except the lower bound on 
the depth of stringer are active. 



(c) When all the constoraints are considered the constraint 
on natuial frequency, panel buckling strength and 
direct stress have become active at the optimijm point 
in all the cases. The skin thickness has not reached 
its lower limit in any of the cases. 

In the case of rectangular ring - rectangular 
stringer combination, side constraint on the stringer 
thickness as well as stringer tuckLing (as column) cons- 
traint are active. For Z - stringer- tubular ring combina- 
tion the thicknesses of rings and stringers have readied 
their respective lower bounds at the optimum point. For 
the combination of Z - stringers and I - rings, the 
lower bound on ring thickness and upper bound on ring 
spacing have become critical at the optimal solution. 

The weight of the stiffened shell has been found 
to have a higher value, when all the constraints are con- 
sidered. This is to be expected as the feasible design 
space becomes narrower with the addition of constraints. 
The stringers are found to contribute more to the total 
weight of the shell than that of rings, with allcons- 
traints. 

5.1.2 Effect of Type of Stiffeners 

As stated earlier, the following stiffener 
combinations are considered for studying the effect of 



type of stiffeners on the optimization problem; 

Case 1 : rectangiolar stringers and rectangular rings, 

Case 2; Z - section stringers and tubular rings, 

Case 3s 2 - section stringers and I - section rings. 

In order to have a meaningful ccanparison of the 
results, the ring cross sectional area and stringer cross 
sectional area are taken to be same in all the three cases 
at the initial design vector. The lower bounds on thick- 
nesses and depths are also taken such that the stiffener 
areas are constrained by the same lower botind in all the 
cases . 

It can be seen that, for all sets of constraints 
considered, I - ring combination (Case 3) gives the minimum 
weight. When the behaviour constraint is only on natural 
frequency, the minimum weight obtained in cases 1, 2 and 3 
are 79 *08, 85.38 and 73-12 Kg. respectively. Thus, the 
minimum weights corresponding to cases 1 and 2 are 8.1 5^ 
and 16.75 % higher than that of case 3} respectively. 

When overall buckling only is considered, the 
minimum weights corresponding to cases 1, 2 and 3 are 
78.3^, 76 . 81 , 76,27 Kg, respectively. Thus cases 1 and 
2 led to 2.72 % and 0.7 1 % higher weights, respectively 
compared to case 3. 



When all the constraints are considered, the 
minimum weights obtained are 189 . 9 ^, 187,61 and 177.77 kg 
in cases 1 , 2 and 3 respectively. This means that the 
minimum weights in cases 1 and 2 are higher by 6.86 % 
and 5*5^/^ respectively compared to case 3 » 

5 . 1.3 Comparison with the Results Reported in Ref. ( 26 ] 

Morrow and Sclmit presented the c^timization 
results of a cylindrical shell \diose data is given in 
Table 5.6, by considering minimization of weight as ob- 
jective with overall buckling and local buckling as 
constraints. In this reference the buckling load has been 
obtained by solving the equilibrium equations. 

The same problem is solved with constraints on 
natural frequency, overall buckling, local buckling, direct 
stress and side constraints. 

The starting design vector is taken to be same 
as that of Ref. TsbJ , and the results of optimization are 
given in Table 5*7. It can be seen that the present 
constraint set resulted in an increase of 31.22^ in the 
minimum weight. A comparison of the present design vector 
with tiiat of Ref. £ 26 J shows that the two design vectors 
(e»‘pecially -the thicknesses of shell) are different at 
■Qie optimal solution. Biis, along with the results of 
Sections 5.1.1 and 5.1.2, shows that the shell thickness 
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tiies to reach its lower hound, value' at the optimum point 
v^en one behaviour constiaint only is considered in the 
problem. Hov/ever, i«iien all the behaviour constraints are 
considered the shell thickness will be in general different 
from its bound value. 

5.2 CONICAL SHELL 

A conical ^ell with rectangular stiffeners, 
having the same overall dimensions as the viking shell L4ol 
is considered for the optimum design, !I5ie optimiEQ design 
of the conical shell is found under three different cons- 
traint sets, !nie first constraint set consists restric- 
tions on natural frequency and design variables. The 
second one contains overall buckling and side constraints 
on design vector. In the third set, constraints are placed 
on natural frequency, overall buckling strength, direct 
stress induced and on the design variables. 

The data considered for numerical confutation is 
shown in Table 5 . 8 . The results of optimization are shown 
in Table 5 . 9 * The computer time taken for the solution of 
one problem ranged from 7 to 12 minutes on lEM 70 M+ 
computer. 

It can be seen that these results ^ow a similar 
behaviour as in the case of cylindrical shells. Some of 
the characteristics are as follows: 
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(i) the minimum weight obtained in the case of third 
constraint set is higher than that of the first 
two cases where only one behaviour constraint 
was considered* 

(ii) the shin thickness approached its lower bound 
value in the first two cases and not in the third 
case, 

(iii) the respective behaviour constraint has become 
active at the optimum point in the first two 
cases, while the constraint on direct stress has 
been found to be critical in the third case. 

In this problem, the ninaber of rings and stringers 
are taken as design variables instead ox stiffener spacings. 

It can be seen that the number of rings reached its minimum 
permissible value in the first case only, -vdiereas the niimber 
of stringers has not reached its lower bound in any of the 
three cases. 

A few optimum solutions of cylindrical and conical 
shell examples are considered for checkin^ the satisfaction of 
Euhn- Tucker conditions and the Lagrange multipliers obtained 
for the solutions are found to be positive. This shows that 
the solutions obtained are atleast local minima. 
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TABLE 5,2 

Nmerical Data for Cylindrical aiell 


Mateilal of the shell, stringers and rings 
Length of the cylinder (a) 

Mean radius of cylinder (R) 

Axial load 

Yield stress of the material 
Factor of safety for the yield stress 
Factor of safety for the buckling load 
Lower bound on natural frequency 
Young's modulus of material (E) 

Poissons ratio (/--') 

Specific weight of the material i ^ ) 
Acceleration due to gravity (g) 

Position of stringer 
Position of rings 


Aluminium 
^3.84 cm 
121 «92 cm 
3.63 X 10^ kg 
2100 kg/cm^ 

1.5 

2.0 

70,0 c/s 

7.35 3c 10^ kg/cm^ 
0.32 

2.713 X lO”^ kg/cm^ 

9.81 nv'sec^ 

Outside 

Outside 


For Z - section web is perpendicular to flange and web to 
flange ratio is 2,5. 

For I - section top and bottom flanges are equal. Web 


to flange ration is 2,0 



TABLE 5.3 


Variable/ § 


obj ective * Bounds 
function | 




I Starting | 
X design “ 
I vector 

J I 


Qr>tlmum design vector. 

^ ‘ witii over-5 -with all 


with f re- 5 
quency 5 
onstraint 5 

i 


all buck- 5 cons- 
ling cons-5 traints 
traint 5 


-t 

0.1270 

0.5080 

0.1275 

0,1270 

0.3033 

w 

■^1 

0.3810 

1,2700 

0.3929 

0,3840 

0 . 4 183 

tp 

0.7620 

1.2700 

1 .^765 

0.7661 

1.5138 

0.3810 

1.2700 

0, If 135 

0.3830 

0.3866 

bp 

0.7620 

1.2700 

0.7798 

2.9726 

1.2997 

1 

2 N-. 381 f 0 

(upper) 

8.1280 

2N-.3590 

24.0280 

19.0250 

d 

76. 6060 
(upper) 

25.5520 

7V. 6760 

75.9460 

13.3860 

f 

•• 

390.7600 

79.0800 

78.^00 

189.9500 

Wi 


257.9500 

64.8000 

64,5200 

153.9600 



100.7600 

12.0900 

6.2100 

16.9100 

Wo 

- 

32.0500 

2.1900 

7.6100 

19.0700 

Frequency ( e/s ) 

78.5807 

70.1584 

- 

70.2400 

Buckling 
strength (Kg) 

1 . 86 x 10 ^ 


5 

7.31x10^ 

8 . 5 x 10 ^ 

Panel 

Ocritical 

1 . 065 x 10 ^ 



1 . 385 x 10 ^ 

■buckling 0 





stressp 

(kg/ cm ) \ induced 

8 . 26 x 10 ^ 



1 . 385 x 10 ^ 

Stringer 5 
buckling 0 critical 

If. 

1 . 20 x 10 



7 . 14 x 10 ^ 

(as CO- 





p 

lumn) 
load kg 

I induced 

1 , 3 ^x 10 ^ 



7*00x10 

Stringer 0 critical 
buckling 5 

33,67x10^ 



2.98x10 

(as plate) 





I 

8.26x10^ 



1.38x10^ 

Direct 

stresso 


8.26x10^ 



1 .38x10^ 

(kg/cm ) 






^7 


TABLE 5, If 


Variable/ 
objective x Bounds 
function ^ 


Optimum design vector 


V VI j. oil j.re- y with over- 0 With all 

vector S quency 0 all buck- | cons- 

il cons- fi ling cons- il traints 
S traint 11 traint 5 


% 

0,1270 

0.5080 


0*1016 

0.25^ 

^1 

0.5080 

1 . 14-30 


0.1016 

0.254-0 

*>2 

1.7272 

3*8100 

1 

2 ^* 38 JfO 

(upper) 

8,1280 

d 

76.6060 

(upper) 

25*5520 

f 


390*7600 



257*9500 

«2 


100 .7600 

’^3 


32.0500 

Frequency c/s 

75-0100 

Buckling 


6 

strength 

(kg) 


4 -, 04 -x 10 

Panel buck-fcritical 

1064- *7 500 

ling 

5 


stress^ 
(kg/cm ) 

1 induced 

825.8100 

\ 

Stringer 

0 critical 

17 . 76 x 10 ^ 

buckling 

5 


(as column^ i^ducpd 
load (Kg) 5 

1 * 34 - x 10 ^ 


Direct 

stressp 825*8100 

(kg/cm"^) 


0.1273 

0.1283 

0.2187 

O.IO16 

0*1052 

0.1017 

1.44-02 

0.5108 

1.2256 

0.1016 

0.1052 

0.1184- 

1.7399 

4- *3099 

5*44-91 

24- *1050 

24-. 3080 

23.0890 

71.044-0 

75*514-0 

9.5760 

85*3600 

76.8100 

18,7.610 

64-. 5600 

65*1600 

111.0500 

18.6700 

6.3200 

16,5000 

2.1300 

5.3300 

60.0600 

70.6800 

- 

71.0000 

- 

7.28x10^ 

7.4-7x10^ 

- 

- 

14-05.0900 

- 

mm 

1399*4^99 

- 

- 

3.46x10^ 

- 

- 

1 .600x10^ 

mm 


1 .394-xlO^ 
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TABLE 5.5 


Variable/ i 

ob j ective I Bounds 

function ^ 

0 


0 _ , 7^ 0 Q-ptimum design vector 

5 Starting 5 vith f re- 0 vith over- d with aiJ 

5 0 quency 5 all buck- 0 cons- 

$ vecior j cons- 0 ling cons- 0 traints 

5 ^ traint S traint S 



0.1270 

O.508O 

0.1270 

0.1273 

0.2235 

*1 

’^1 

0.1016 

0,2540 

0.1016 

0.1019 

0.1054 

1 . 53 (A- 

3.4290 

1 .6203 

1.5578 

1.5692 

4 

0.1016 

0.2540 

0.1029 

0.1039 

0.1359 

^2 

1.7272 

3.8100 

1.7871 

4.2172 

4.6703 

1 

24.3840 

(upper) 

8.1280 

24.1810 

23.8760 

, 23.4950 

d 

76.6060 

(upper) 

25.5520 

74.3460 

72.5930 

9.7540 

f 

- 

390 .7600 

73.1200 

76.2700 

177 .7700 

^2 

- 

257.9500 

64.5200 

64.6100 

113.4900 

- 

100,7600 

6 .4900 

6 .3000 

6 .6800 

4 

- 

32.0500 

2.1100 

5.3600 

57.6000 

Frequency c/ s 

121.5400 

70.6700 

- 

74.5200 

Buckling 
strength (kg) 

4 . 45 x 10 ^ 

- 

7 . 77 x 10 ^ 

4 . 4 lx 10 ^ 

Panel 5 critical 

1064.7500 



I 4 l 4 .7700 

buckling y 





1400.0000 

stress^ 0 induced 

825.79 

- 

- 

(kg/cm ) 5 






Stringer 

fcritical 

177210.13 

- 

-- 

23591 .7600 

bucko ing 






(as column^ 
load (kg) 5 

1337.76 

- 

- 

1554.4800 

Direct 
stressp 
(kg/cm ) 


825.7900 

- 


1400.0000 


TABLE 5.6 

Nime ideal Data for the Example From Bef, T 26 I 


Length of cylinder (a) 

Mean radius of cylinder (R) 

Yield stress of the material 
Poissons latio 

Lower bound on natural frequency 

Young’s modulus of material (E) 

Axial load per unit length of 
circumference 

Specific \ieight of material ( ^ ) 

Factor of safety for direct stress 
Factor of safety for buckling strength 
Placing of stiffeners 
Type of stiffeners 


96.5 cm 

^.2 cm 

3.5 X 10^ kg/cip^ 
0.33 

100 c/s 

7,35 X 10 "^ kg/cm 
lk-3.2 kg/ cm 

0.00278 kg/cm^ 

1.5 
2.0 

Outside 

Rectangixlar 
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TABLE 5.7 


Variable/ x 
oh 3 ective i 
function ^ 

Bounds ■ . 

T 

5 

0 

5 

starting 

design 

vector 

Q Ontimum design vector 

5 with present 5 As given 

5 constraint 0 in _ 

5 set 5 Ref. (26} 


O.0516* 


0.0711 

0.0442 

0.02360 


0.025+* 


0.1270 

0.0254 

O.0516O 


0.1270 


0.2540 

0.1748 

4.64800 

^2 

0.0254* 


0.1270 

0.0488 

0.03810 

^2 

0.1270 


0.2540 

0.5286 

O.36O70 

1 

12.7000 

(upper) 


3.8100 

6.0096 

6.98500 

d 

5.0800 

(upper) 


0.6350 

1 .5037 

O.67O6O 

f 

- 


5.3600 

2.5589 

1 .95000 


~ 


2.9300 

1 .8222 

0.96500 


- 


0.3500 

0.03157 

0.97600 

W3 

— 


2,0860 

0.70520 

0.00900 


* These hounds in the original Eef [26j are 0,0 



51 


TABLE 5.8 

Nuanerical Data for Conical ShPill 


Material of shell, stringers and rings 
Mean radius of cone at the small end (R^) 
Mean radius of cone at the base (R2) 

Cone apex angle (90 -c^) 

Axial load 

Yield stress of the material 
Factor of safety for yield stress 
Factor of safety for buckling 
Lower bound on frequejicy 
Youngs modulus of material (E) 

Poissons ratio (/-<-) 

Specific weight of material ( ^' ) 
Position of stringer 
Position of rings 


Magnesim 
81 .28 cm 
175*26 cm 

70.0° 

If. 536 z 10^ kg 

1120 kg/cm^ 

1.5 

2.0 

75.0 c/s 

6 2 

if ,06 X 10 ^ kg/cm 
0.33 

1.8279 X I0"^kg/cm“' 

Outside 

Outside 
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TABLE 5 

.9 



Variable/ 

Objective 

function 

1 Lower 

1 bounds 

0 

5 

1 Starting 5 
0 design J 
0 vector 5 
0 0 
D 0 

Optimum 
With fre-^ 
quency 0 
cons- 5 
traint I 

design vector 

all buck- 5 cons- 
ling cons-0 traints 
traint 5 

% (cm) 

0.1270 

0. 508O 

0.1280 

0.1331 

0.1839 

t^ (cm) 

0. 3810 

1. 2700 

0.3924 

0.4013 

0. 4648 

b^ (cm) 

0.7620 

1. 2700 

0.7628 

0.7874 

0.7945 

±2 (cm) 

0.3810 

1. 2700 

0.5820 

C.5OO4 

0.6761 

b^ (cm) 

0.7620 

1. 2700 

0.7884 

5.3421 

5.8473 


10 

30 

12 

23 

28 


10 

30 

18 

17 

21 

f (kg) 

- 

155. 2100 

25.1730 

35.5705 

57.5775 

(kg) 

- 

74.9600 

18.8800 

19.6431 

27.1194 

¥2 (kg) 

~ 

71.3900 

5.3000 

10.7224 

15.2555 

¥^ (kg) 


8.8600 

0.9930 

5.2050 

15.2026 

Preq.uency c/s 

104.4800 

79.9100 

- 

199.1500 

Buckling 
strength (kg) 

3.17x10^ 

- 

1.95x10^ 

11.74x10^ 

Direct stress 
(kg/cm^ ) 

6.13x10^ 

- 

- 

10.662:10^ 



CHAPTER 6 


CONCLUSIONS AND EECO^i^IENrATIONS 
6.1 CONCLUSIONS 

Hie following conclusions can be drawn from the 
results obtained for the optimization of cylindrical and 
conical shells; 

i) When the natural frequency or overall buckling 
strength is the only behaviour constraint con- 
sidered, it is becoming active at the optimum 
point. 

ii) The objective function (wei^t of the ^ell) is 
found to be 'very sensitive to the design variable, 
t^_^. The side constraint on t^ is becoming 
active when only one beliaviour constraint and the 
side constraints are considered. 

iii) In most of tlie cases the side constraints on thick- 
nesses of ring and stringer are becoming active. 

iv) Cut of the three combinations of stiffeners con- 
sidered, the Z - section stringers and I - sec- 
tion rings combination is found to give minimum 
weight . 
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r) In the case of conical shell, the numher of rings 
is coming to the lower bound under any constraint 
set, 

vl) The optimum weights of the cylindrical as well as 
conical shells are found to increase as tiie number 
of constraints considered are increased, 

vii) The stringers are seen to contribute considerably 
to tile total weight of the shell, when all the 
constraints are considered. 

6.2 KSCOM'ISNDATIONS FOR FUTUK3 VJOEK 

i) The procedure adopted in forming the natural fre- 
quency and overall buckling constraint is an 
approximate one. It can be made more accurate by 
searching for the minimum frequency and buckling 
strength in the (m, n) space at each design vector 
in the optimization routine, 

ii) The influence of inside stiffeners or combination 
of inside and outside stiffeners (rings inside - 
stringers outside or stringers inside - rings 
outside) on the optimum design can also be stu- 
died. 

iii) In order to study the effect of type of stiffeners 
on the design of conical shells, a number of 



combinations of stiffeners are to be considered, 
in addition to the rectangiG.ar stiffeners con- 
sidered in this work. 

The design optimization can be carried witii other 
criteria like maximization of minimtm natural 
frequency and maximization of the difference bet- 
ween the lowest two natural frequencies. 

The mixed-integer nonlinear programming techni- 
ques can be used to solve the problan with 
integer variables (namely, number of rings and 
number of stringers) more efficiently. 
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APPENDIX A 


Expressions for of Eqs. (3*12) and (3*17) 

= D (IN8) n + D^ - n^ + . 

w cos'^ r cosc; 2 

D^ IT coso( j_ (IN8) + 2 kjjj (IN lif) + (IN 4)1 

+ (IN 7) -y ?- ■ (IN 6) 


cos o< cos M 

rt- - 2 , o (1 - P ) /Tk^ sin\- 

+ 3 (IN 6) n + 2 SHTT 

cos ==s 


(IN 2) 


+ T. 


o 2 

k'^ sin c< 

m D 


cos ->' 


(IN 1) (IN 16) 


- \ n (IN 13) n - D^ /T (IN 13) n - D^N '2 (IN 13) 

- kj^ (IN 9)] n + \n [(IN 13) + kjjj (IN 3)] n 


Sp '2> " 


P . (IH 8) n - D^ . C . s . 3-ng . ^ 

W cos ty. r cos:^' 


2.2^ 

+ Sr /tUiN 7) (2 - ^ S - m .. ri -) + (IN 13)1 n 

cos ^ 


+ K 


+ K 


/f sini:>; r, 


v1 


W cos Of 


^ I (IN 6 ) (2 i^) + kjj^ (IN 12 ) I n 

~ cos OK 

i l mng _ r (IN 6) (2 - -I- ) 


r coso*. L 


cos ^ 
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+ If (IN 12)1 n + K 

m J W w cosoe 


2 (IN 6) 


+ 4 ]k^ (IN 12) - (IN 8) I n 
2 /r k^ sino; 


(1 -7%) 
w 


K 


w 


COSO# 


i (IN 12) + kjj^ (IN 2)] n 


2 

k sino{ 

- -”^cosoc (IN (IN 16) n 


r>^ n cosof (IN 2) + D^/r coso? (IN 2) 


+ rr COSO*, f 3 (IN 2) - 4 kjj^ (IN 14) + k^ (in 10)] 

+ Dg cos:;i I 4 (IN 1) - 4 kjj^ (IN 13) + k^ (IN 9) ] (IN 16) 


+ /^ 2 /r cosc^ 12 (IN 2) - kjjj (IN 14) 


(1 - Hr) rf 5 

^ \ <™ 2> “ 


N" sinc< (IN 13) + H sino\ (IN 13) 

+ /^ /r sin=^ [2 (IN 13) - kjjj (IN 9)] 

- /r coso\ [(IN 12) (2 ) + kj^ (IN 2)1 

COS ^ ^ 

r 

- Sg coso^ 4 (IN 11) + kj^ [ 8 (IN 1) - 2 (IN 7)] 
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H 


6 = 


D J l . sin ^ /I 3221 ^ 0 ^ ^ 

w cos + -^r cos oi 


- S^ 2 /r since [(IN 7 ) (2 - -i£- ) + ^ (IN 13)i 

cos C! “ 


r 2 

+ n- cosc^ I (IN 6) (2 % — 


+ 2 (IK 12) (2 - ^ii-) (IN 2)1 

cos 


2 

+ cosc( |~(IN 6) (2 + 2 k (IN 12) 

onQ rk ^ 


2 

(2 (IN 2) J + /T cos jS< | ^ (IN 6) 


2 

cos CA 


+ 


16 kj^ (IN 12) - if k^ [ (IN 8) - If (IN 2)1 


- 8 k^ (IN llf) + k^ (IN 10) I + Kg cosoi If (IN 5) 

I 

+ 16 kj^ (IN 11) - if k^ f (IN 7) - if (IN 1)] - 8 k^ 

+ (IN 9) I (IN 15) + 2/rcos.o(j r2(IN 6) 

J ’! ^ 

s. 

2 

+ if kjj^ (IN 12) - k^ (IN 8)j (2 % ) 

2 k^ (IN 12) + ^ k^ (IN 2) - l5^ (IN lif ) I 

[(IN 6) + 2 (IN 12) 


(1 -/i:j 

iL. N 

2 -w coso^ 


k^ (IN 2) I + T 


/r 


r cos o 


if (IN 6) 


+ if kjjj (IN 12) + k^ (IN 2) 


+ T 


s cos oi 


(IN 1) 


(IN 16) n' 


13) 
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^7 cos<^ (IN 10) 

+ nig cosd. (IN 9) (IN 16) 

^8 ~ oosc\ (IN if) + m^/r coso\ (IN if) 

+ nig cos^ (IN 3) (IN 15) 

/r cos^( (IN 10 ) + m^rj coso^ (IN 10 ) 

+ m„ coscH. (IN 9) (IN 15) 

w 

^10 "" ^^1 r 3) + ^ (IN 7) + ^ kj 2 (IN 13) j cot o{ 

■where 

(IN 1) = ^ (S^ - S^) + (1/lf 1?^) (s^ - S^) 

(IN 2) = ^ (S^ - S\) + (3/8 k^) (S^ - S^) 

(IN 3 ) = ~ (S^ - S^) + (1/2 k^) (S^ - S^) - ( 3 A k^) (S 2 - S^) 

(IN h) = ^ (s| - S^) + (5/8 k^) (S^ - S^) - (15/8 k^)(s| - S^) 
(IN 5) = ^ (S 2 - s^) 

(IN 6) = ^ (S| - S^) 

(IN 7) = ^ (S| - S^) - (1/if (S 2 - S^) 

(IN 8) = ^ (S^ - S^) - (3/8 k^) (s| - -S^) 

(IN 9) = ^ (S^ - S^) - (1/2 k^) (s| - S^) + ( 3 A “ ^1^ 

(IN 10) = T2 ^4 “ 4^ “ ^4 ’ 4^ “ 4^ 
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(IN 11) = - (1/f k^) (s^ - S^) 

(IN 12) = - (1A l5^) (S^ - S^) 

(IN 13) = - (1A (s| - Sp + (3/8 k^) (S 2 - S^) 

(IN A) = (1A k^) (8^ - s^J-) + (3A k^) (s| - Sp 

(IN 15 ) = /T and 

(IN 16) = rr 



